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BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates generally to multiple access spread spectrum 
radio receivers, and more particularly to receivers with enhanced ability to acquire 
and track a relatively weak signal in the presence of a comparatively stronger signal. 

The difference in signal strength often can be attributed to the relative 
distance of the signal source and the receiver, and thus the difficulty in tracking the 
weaker signal in the presence of a closer, stronger signal is often referred to as the 
near-far problem of spread spectrum multiple access. This problem can also occur 
when one signal source is obscured from the receiver while another signal source 
has a direct line of site. An example of this would be operating a receiver inside a 
building, perhaps near a window oi- a door, thereby receiving some signals at normal 
signal strength while others are attenuated by the building structure. 

Description of the Prior Art 

The Global Positioning System (GPS) is a radio navigation system operated 
by the United States Air Force for the dual purpose of providing accurate global 
positioning information to military as well as civilian users. To this end GPS provides 
two services: the Precise Positioning Service (PPS) which is available primarily to 
the US arrned forces and requires the use of receivers equipped with the proper 
PPS equipment, and the Standard Positioning Service (SPS) which is less accurate 
than PPS but is available to all users whether or hot they have access to PPS 
equipment. The U.S. Department of Defense has the capability to degrade the 
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accuracy of the SPS through what is known as the Selective Availability (S/A) 
algorithm, and has taken an official position that all such S/A induced errors will be 
limited to a 100 meter horizontal position error range (2d-RMS). In contrast, PPS is 
accurate to within 22 meters. 

GPS is essentially comprised of at least 24 satellites in orbit around the Earth 
at an altitude of approximately 20,000 Km in one of six orbits; Each orbit is occupied 
by at least four satellites. Each GPS satellite broadcasts a unique radio ranging 
signal that can be received by properly equipped GPS receivers. The signal 
contains information identifying the particular transmitting satellite and navigation 
data such as time and satellite position. On a fundamental level, all GPS receivers 
operate by tracking the ranging signals of multiple GPS satellites and determining 
the user's position in terms of latitude, longitude, and altitude or another equivalent 
spatial coordinate system. 

The ranging signal broadcast by each satellite is comprised of two signals: 
the primary Link 1 (L1 ) signal broadcast at a carrier frequency of 1 575.42 MHz and 
the secondary Link 2 (L2) signal broadcast at a carrier frequency of 1227.6 MHz. 
Both LI and L2 carrier signals are spread spectrum signals modulated by digital 
signals, or codes, that "spread" the spectrum of each carrier signal over a specific 
bandwidth. The LI signal is modulated by three bi-phase (i.e. ±1 ) digital signals: the 
Clear or Coarse Acquisition (C/A) Code which is a short Pseudo-Random Noise 
(PRN) code broadcast at a bit (or chip, which refers to each.pulse of the noise code) 
rate of 1,023 MHz (and thus spreads the L1 carrier signal over a 1.023 MHz 
bandwidth by essentially breaking each bit in the original signal into 1023 separate 
bits, or chips, in what is known as direct sequence spread spectrum) and which 
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therefore repeats every 1 millisecond; the Precise (P) code which is a much longer 
PRN code that repeats every week and is broadcast at ten times the chip rate of the 
C/A code (10.23 MHz); and a 50 Hz navigation data code (D). The C/A code is 
always broadcast in the clear (or unencrypted) whereas the P code is encrypted by 
an encrypting (E) code to form what is known as the Y code. The low data rate 
navigational code D comprises orbital parameters and clock correction information 
for the satellite modified by S/A. 

Currently the SPS is predicated solely upon the L1 signal but in the future the 
SPS signal will be available on both L1 and L2. The current L1 signal contains an 
in-phase component modulated by P©E©D (where ® denotes the logical XOR 
function) and a quadrature component modulated by C/A©D, and can be 
represented for each satellite i as 

Siu(t) = -yJlPlU X ei{t)pi{t)di{t)cos[o)L\t + (pL\] 
+2^ Plu X Ci(t)di{t)sin[^ci)L\t ^- ^ii] 
where A represents the signal power, © the carrier frequency, and cp a small phase 
noise and oscillator drift (i.e. clock error) component. 

The broadcast satellite navigation data message D and algorithms to process 
it are defined in the publicly available U.S. government specification ICD-GPS-200. 
The satellite position portion of D is actually a prediction that is computed using 
ranging measurements of the GPS satellites taken at five monitoring stations 
distributed around the Earth. Periodically, typically daily, the OPS control segment 
uploads each satellite with its predicted navigation data and an estimated correction 
to its on-board atomic clock. 
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The satellite navigation data includes the GPS almanac which is used to 
predict the position and velocity of each GPS satellites for many weeks into the 
future. A typical GPS receiver uses the almanac data, the algorithms defined in ICD- 
GPS-200 and standard linear equation solving techniques to compute the position 
and velocity of each GPS satellite and to predict the expected range (PRN code 
phase) and Doppler frequency at which the receiver will find the satellite's signal. 

Because all satellites broadcast at the same carrier frequency, each of the 
satellite ranging signals must be able to share this frequency with a minimum of 
interference from the other signals. This is accomplished by carefully selecting the 
PRN codes to have a sharp (1-chip wide) autocorrelation peak to enable code- 
synchronization and achieve equal spreading over the whole frequency band, and 
further have low crosscorrelation values, in a method known as Code Division 
Multiple Access (CDMA). The C/A PRN codes are unique to each satellite and are 
taken from a family of codes known as Gold codes. The GPS C/A codes are formed 
as the product (or modulo-2 sum) of two maximal binary code sequences (G1 and 
G2)each 1023 bits long . The 1023 members of this Gold code family are generated 
by shifting the starting state of the G2 register with respect to G1 . Thirty-two out of 
the 1023 possible Gold codes were selected for the GPS satellites based upon two 
criteria: the number of ones and zeros in the code must differ by exactly one (i.e. the 
codes are balanced), and the crosscorrelation between any two of the C/A codes 
is no more .than 65/1023 or -23.9 dB (normalized to the autocorrelation peak of 
unity). This crosscorrelation immunity is called the Gold bound, and represents the 
maximunh interference between equal strength C/A code signals with identical 
frequencies. This PRN signal design enables satisfactory CDMA operation of the 
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GPS system, i.e. as many as 32 satellites sharing the same broadcast band, 
provided that the received powers of the GPS signals are not larger than the Gold 
bound, which is typically the case. 

The Gold code bound is applicable for signals with identical carrier 
frequencies. However, due to Doppler frequency shifts caused by motion of the 
satellites in their orbits and movement of the receiver, the received frequency of the 
GPS satellite signals is typically shifted by up to ± 5 KHz from the nominal 1 575.42 
MHz LI carrier frequency. Relative to any single satellite, the frequency of other 
satellites may differ by as much as ± 9 KHz. 

The strong/weak crosscorrelation problem is worse if the signals are Doppler 
shifted. As mentioned previously, the C/A code's Gold code family is generated by 
forming the mod-2 sum of a selected pair of maximal binary code sequences (G1 
and G2) for all 1023 possible time shifts between the two sequences. The 
crosscorrelation (which is the multiplication of two signals) for binary codes is 
equivalent to mod-2 addition of the codes because multiplication of ±1 values has 
a one-to-one correspondence with mod-2 addition of binary 0,1 values. Therefore, 
the crosscorrelation of two Doppler shifted members of the Gold code family 
reduces to mod-2 addition of each maximal sequence with itself, followed by 
another mod-2 addition. The shift-and-add property of a maximal sequence means 
that the mod-2 sum of a maxihial sequence with a shift of the same maximal 
sequence yields yet another shift of the same maximal sequence. Therefore, 
crosscorrelation of two Doppler shifted members of the Gold code family yields 
another member of the same Gold code family. It has been found that these 
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generated Gold codes are not mennbers of the C/A family and may have 
crosscorrelations that exceed the C/A code design limit. 

No closed form analysis of the crosscorrelation interference of Doppler 
shifted C/A codes with relatively different carrier frequencies is known. Instead, 
simulations are used to analyze the effects of Doppler shifts on the crosscorrelation 
of C/A codes. The simulations either generate the two desired frequency offset 
codes and compute the crosscorrelation directly or generate the Fourier transform 
of each code, adjusted for frequency offset, and compute the crosscorrelation of the 
transforms. It has been found that for a ±9 KHz Doppler range the worst case 
crosscorrelation for the GPS C/A codes is -20.9 dB. This worst case scenario occurs 
when the relative Doppler shift between the two satellite signals in an integer 
multiple of 1 KHz. 

While Doppler offsets increase the level of strong/weak signal 
crosscorrelation when the frequency difference is an integer multiple of 1 KHz, 
frequency attenuation decreases the crosscorrelation effects when the Doppler shift 
is not a multiple of 1 KHz. The GPS receivers integrate (sum) the in-phase and 
quadrature (I, Q) measurements for some length of time before they are used for 
signal detection or signal tracking. If the integrated signal contains a frequency 
error, then the accumulation decreases the signaPs apparent strength by the well 
known si n(x)/x function, where x is half the amount of phase rotation in radians that 
occurs over the integration period (note that the limit of sin(x)/x is 1 as x approaches 
0). Thus, if the Doppler difference between the replica weak signal and the 
interfering strong signal is 500 Hz and the I. Q integration time is 1 ms, then x is 
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equal to 7i/2 radians, sin(x)/x is equal to 2ln and the interference is attenuated by 
approximately 4 dB. 

Consequently, a strong/weak signal crosscorrelation problem may occur if the 
strength of one satellite approaches being 20.9 dB stronger than the strength of the 
second satellite. Under this condition, the acquisition search may detect the 
crosscorrelation spectral line from the strong satellite instead of the autocorrelation 
spectral line from the weak satellite. 

The GPS system was designed with the assumption that receivers would be 
operated out-of-doors with direct lines of sight to all satellites. In this case the C/A 
code provides adequate protection against strong/weak signal crosscorrelation. 
However, once a receiver moves indoors or under a canopy of trees, some of the 
signals can become significantly attenuated while the others continue to be received 
at normal signal strength. In such circumstances the operational significance of the 
crosscorrelation peaks of the Gold codes is to cause difficulty in being able to 
discriminate between a weak GPS signal and the crosscorrelation of a relatively 
stronger GPS signal. An incorrect discrimination may cause large errors in the 
latitude, longitude and altitude computed by the GPS receiver. 

An SPS-equipped GPS receiver will receive at any given time the L1 ranging 
signals from as many as twelve satellites, all multiplexed on the same carrier 
frequency, each modulated by its own C/A PRN Gold code. From this compound 
carrier signal the receiver must be able to identify and extract the individual 
satellites' signals and then process each of these signals to recover the information 
contained therein. Each of these satellites has the potential of interfering with every 
other satellite signal. In a worst case, when the signals from a single weak satellite 
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and a plurality of strong satellites are received simultaneously, the weak satellite 
signal may have significant crosscorrelatlon interference from each strong satellite 
signal. 

When a GPS receiver is first powered on, it has at best only an approximate 
knowledge of its position, its local oscillator offset (which will appear as a Doppler 
frequency offset that is common to all satellites) and the correct time. Therefore, the 
receiver must perform a systematic search through a large portion of all possible 
C/A code phases and all possible Doppler offsets to locate the satellite signals. 
During the search the strong/weak crosscorrelatlon from any relatively strong 
satellite may cause the receiver to mistake a crosscorrelatlon spectral line from the 
strong satellite as a signal from a weak satellite. 

After the receiver has started it can predict the C/A code phase and Doppler 
offset of all the satellites using the almanac data and the algorithms of ICD-GPS- 
200, ^at which time it only needs to search a relatively smaller range of C/A code 
phases and Doppler frequency offsets for the desired satellite signal. Nonetheless, 
the strong/weak crosscorrelatlon problem remains when the crosscorrelatlon peak 
from a relatively strong satellite occurs within the search range of a relatively weak 
satellite. 

A typical GPS receiver consists of an antenna to receive the carrier signal 
while rejecting multipath and. optionally, interference signals; a preamplifier 
comprising a bandpass filter to filter out potential high-leveL interfering signals in 
adjacent frequency bands, and a low noise amplifier (LNA) to amplify the carrier 
signal; a reference oscillator to provide time and frequency reference for the 
receiver; a frequency synthesizer driven by the oscillator; a downconverter to 
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convert the filtered carrier signal to an intermediate frequency (IF); an IF section to 
provide further filtering of out-of-band noise and interference, amplification of the 
signal to a workable signal-processing level, and optionally down conversion of the 
IF signal to a baseband signal; and an analog to digital converter (ADC) to sample 
and quantize the signal into in-phase (I) and quadrature (Q) components. The ADC 
may sample either the IF or the baseband signal, depending upon the receiver 
design. 

The digitized I. Q signal is next fed into one to twelve or more tracking 
channels. There it is correlated with a C/A PRN code replica that may either be 
generated internally as needed with shift registers or stored as a complete set of 
precomputed C/A code chips in memory. The correlation of the replica and received 
signal is accomplished by mixing (multiplying) the two signals and integrating 
(summing) the power of the in-phase and quadrature components of the resultant 
signal. Typically, the phase of the carrier and code of the replica signal are aligned 
with the received signal by Phase Locked Loops (PLLs), Costas Phase Detection 
Loops and/or Delay Locked Loops (DLLs). The PLL and Costas loops maintain 
phase agreement of the received and replica signal by driving the quadrature 
component to zero while maximizing the in-phase component. The DLL maintains 
C/A code alignment by balancing the correlation power measured at two or more 
code offsets such as early and late or early and punctual. Each recovered spread- 
spectrum LI signal is then fed to the signal processing portion of the receiver where 
it is demodulated to recover the signal carrier and the C/A and D codes. These in 
turn are supplied to a navigation data processor where the position of each satellite 
being tracked is computed from the D code and various error corrections are 
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performed. Sources of error include ionospheric and tropospheric delays, the 
Doppler effect, satellite and receiver clock errors, equipment delays, noise, and 
multipath errors due to a signal, being reflected and thus received multiple, but 
slightly delayed, times. 

The maximum C/Nq (signal to noise ratio in a 1 Hz bandwidth), received by 
GPS receivers near the surface of the Earth is approximately 55 dB-Hz. allowing for 
additive multipath interference. In contrast, state of the art GPS tracking algorithms 
can acquire and track GPS signals with C/Nq as low as 24 dB. and future advances 
promise to lower this threshold even further. Thus, the range of useable GPS signal 
power is 35 dB or more. Assuming a worst-case strong/weak crosscorrelation C/A 
code spectral line of -20.9 dB, a method is needed to increase the discrimination of 
the C/A code by at least 10 dB-Hz. 

The prior art has developed a general approach for predicting the 
crosscorrelation of two Doppler shifted PRN code sequences if the code timings, 
carrier phases, and signal amplitudes are known. The solution may be summarized 
as optimum maximum-likelihood demodulation of the unknown data bits by a 
computationally intensive Viterbi algorithm. In a practical sense, this optimum 
demodulation can be viewed as equivalent to doing strong signal cancellation with 
enough delay introduced to estimate the unknown data bits of the strong signals 
with a low error rate. This general solution assumes an ideal channel, but a practical 
solution to the near-far problem of CDMA must cope also with multipath propagation 
effects. 

Hence, in light of the above, those skilled in the art have recognized a need 
to increase the strong/weak signal discrimination within CDMA coded spread 
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spectrum signals. It has also the been recognized by those skilled in the art that it 
would be of value to develop such a method that will be compatible with the SPS of 
GPS. The present invention satisfies these needs as well as others. 

SUMMARY OF THE INVENTION 

What has been needed and heretofore unavailable is a method for removing 
the effects of a strong code spread signal on a weaker code spread signal, the so 
called near/far or strong/weak problem of CDMA, which can be implemented in 
existing systems without exceeding system throughput limitations. 

The method of the invention allows for a post correlation removal of strong 
signal effects on a weaker signal, and can be implemented in almost any 
multichannel receiver with only a modest addition to the overall throughput loading. 
The resulting corrected weak signal extends the operation of CDMA receivers into 
traditionally difficult areas such as in and around buildings or under a forest canopy. 

In general, the method consists of tracking one or more strong signals in a 
multi-channel CDMA receiver such as a GPS receiver. Using information about 
available signal sources, the receiver may classify any signal sources believed to 
be present but not currently being tracked as weak signals. These weak signals may 
be tracked-by removing the crosscorrelation effects of all the-strong signals on the 
weak signals. This is done by setting a channel of the multi channel receiver to the 
predicted frequency and code phase of each weak signal. The measurement from 
this channel will contain the crosscorrelation of any strong signals with the desired 
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weak signal. The crosscorrelation can be calculated by crosscorreiating the code 
sequences of the strong and weak signal channels. Because the strong signal Is 
being tracked, its amplitude and phase are known. Finally, as discussed previously, 
the crosscorrelation has nnaxinnum peaks when the relative Doppler between the 
signals is an integral multiple of 1000 Hz. By scaling each tracked strong signal by 
the attenuation caused by the difference in frequency between the strong and weak 
signal and multiplying by the calculated crosscorrelation, the effect of the strong 
signal on the weak signal can be estimated and thus removed. To allow for both 
carrier and code tracking of the weak signal with the PLL, DLL and Costas loops, 
the process must be repeated for at least two reference code offsets such as early 
and late or early and punctual. 

Signal detection of a weak signal can be handled in either of two ways. The 
simplest method entails only performing signal detection when the delta frequency 
(difference in actual received frequencies) between the strong and weak signal 
provides sufficient attenuation of the crosscorrelation of the strong signal with the 
weak signal. The more complete, but slower and more complicated, method is to 
search over the appropriate range of Doppler frequencies and possible code offsets, 
using the method of removing the crosscorrelated strong signal for all possible 
Doppler and code offsets, 

Otherfeatures and advantages of the present invention will become apparent 
from the following detailed description, taken in conjunction with the accompanying 
drawings, which illustrate, by way of example, the principles of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a flowchart depicting the main steps of the method of the 
present invention. 



DETAILED DESCRIPTION 



The strong/weak or near/far signal isolation provided by a spread spectrum, 
PRN code family such as used in CDMA is dependent upon the crosscorrelation 
between the various code members of the family. In the case of GPS, the isolation 
of two signals at the same frequency (or multiples of the code repetition rate, in this 
case 1 KHz) is about 21 to 23 dB. If the relative strengths of two signals differ by 
more than this limit, the weaker signal cannot be discriminated using only the 
spreading code. A method of removing the effects of the stronger signal must be 
applied if the weaker signal is to be tracked. 

As noted above, the crosscorrelation effect is at its maximum when the 
relative Doppler frequency offset between the relatively strong and weak signals is 
an integer multiple of 1 KHz in the case of C/A GPS signals. 

A general solution to the problem of tracking a weak signal spread spectrum 
signal in the presence of a stronger spread spectrum signal has been developed. 
It is based on the premise that all aspects of the strong signal's interference can 
either be measured or calculated in order to remove it from the weaker signal. The 
solution can be implemented in any multi-channel receiver having the ability to 
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control a channel's frequency and phase as well as selecting the desired spreading 
code and setting that code's phase position. The receiver typically employs two 
channels, one to track the weak signal and one to track the interfering strong signal. 
However, the channel that is used to track the strong signal is not required if the 
characteristics such as power, code phase and frequency of the strong signal can 
be obtained or accurately estimated by alternate means. 

As shown generally in Figure 1 , the strong signal is first acquired 10, such as 
by tracking in a first channel of the receiver. The channel provides a measurement 
of the signal strength of the strong signal along with the phase of the carrier signal 
and the spreading code. Additional channels may be used to track additional strong 
signals (not shown in the flowchart). 

The code phase of the spreading code of the weaker signal, along with its 
received frequency and signal phase, are predicted 20 based on data message D 
by methods known in the art. A second channel in the receiver is dedicated to 
receiving the compound carrier signal and tracking 30 the predicted weak signal 
component. 

The second receiver channel correlates the incoming signal with the second 
code at the predicted frequency and signal phase. The resulting in-phase and 
quadrature (I, Q) measurements contain both the weak signal and the strong signal, 
each spread by their unique code. Correlation by multiplication of the replica code 
for the second signal, Code2R, with the incoming signal yields the product 
Code2Rx(weak2xCode2 + StrongX^CodeX + ...) where weak2 is the power of 
weak signal 2, Code2 is the actual code for satellite 2 broadcasting the weak signal 
2. StrongX (X = 1 , 3, 4, ...) is the power of strong signal X, and CodeX is the actual 
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code for satellite X contained in the signal. The product Code2RxCode2 is the 
autocorrelation of the received code 2 and the replica code 2. The autocorrelation 
function has a value of 1 if the replica code is aligned with the received code. This 
crosscorrelation of replica code 2 with code X (Code2RxCodeX) is next computed 
40 to be removed from the compound signal. 

Code1 and Code2 are both members of a PRN code family, and their 
autocorrelation and crosscorrelation properties are known. It is therefore possible 
to calculate the crosscorrelation of the two codes at their respective phases by 
simply multiplying each bit of Code1 by the corresponding (in time) bit of Code 2 to 
produce their crosscorrelation value. Since there may be a relative Doppler 
frequency offset between the two codes, the phase of the codes will precess past 
one another over time and create a new crosscorrelation function. For the GPS 
system the greatest delta code Doppler typically encountered is about ± 9 KHz 
which is equivalent to 6 code chips per second (1 540 carrier cycles per code chip), 
and thus the maximum recalculation rate of the crosscorrelation value is roughly 6 
times per second. 

As previously indicated, maximum crosscorrelation occurs at a frequency 
offset of zero with peaks occurring at intervals of 1000 Hz. There is an attenuation 
of the crosscorrelation as the frequency offset moves away from zero. This 
attenuation follows the well known sin(x)/x curve. If 10ms measurements are used 
for tracking- or acquisition, the attenuation factor would be. equal to sin(Afreq x 
7c/1 OOHz)/(Afreq x tt/IOOHz). This produces an attenuation of-10dB at abouta 75Hz 
delta frequency. Other local peaks in the sin(x)/x cun/e (i.e. locally minimum 
attenuation) occur at 150 and 250 Hz with attenuations of -13.5 and -18 dB, 
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respectively. This implies that for a desired strong signal suppression of 1 0 dB, only 
the first lobe of the sin(x)/x function need be considered; however, should additional 
suppression be desired, the entire curve may be considered. 

The next step entails computing 50. for each strong signal, the product of the 
strong signal amplitude and the calculated frequency and time domain (code phase) 
crosscorrelation. The weak signal is finally extracted by subtracting 60 this product 
from the compound signal. The weak signal thus extract is subsequently processed 
in the receiver circuitry as known in the art. 

The in-phase and quadrature amplitude (I, Q) of each strong signal is 
obtained by measurement in each strong signal's own individual receiver channel 
or by estimation through independent means. Because the strong signal is being 
actively tracked by the receiver's phase lock loops, the phase of the strong signal 
is presumed to be near zero radians and thus nearly all the signal power is in the 
in-phase portion. 

A signal comprising a strong signal S1 modulated with a first code Code1 
summed with a weak signal w2 modulated with a second code Code2 produces (S1 
>^ code1 + w2 X Code2). The sum of the two signals is correlated with a replica of 
the second code Code2R to produce i:{Code2R x (S1 x Code1 + w2 x Code2)}, 
where the sum Z includes all chips of the PRN code used to modulate the weak 
signal w2. The auto-correlation of a code with itself is 1 so the preceding equation 
can be rewritten as S{S1 xCodel xCode2 + w2}. We can see that in order to obtain 
w2 we must remove S1 xCodel xCode2. Since we know Code1 and Code 2, we can 
easily calculate their crosscorrelation. This leaves us to estimate the value of S1 
which can be done by independently tracking the strong signal on a separate 
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channel, or by any other convenient means. This computed value of 
S1 xCodel xCode2 would be sufficient if the strong signal S1 and the weak signal 
w2 where at the same frequency. The two signals are received at different 
frequencies, however, due to the Doppler effect as well as other factors enumerated 
previously. We know that strength of the crosscorrelation varies with the difference 
between these frequencies in a sin(x)/x relationship. We must therefore calculate 
an attenuation factor based on the difference in frequency between the strong and 
the weak signal and apply it to the computed crosscorrelation. Furthermore, if more 
than one strong signal is present, an attenuation factor must be computed for each 
strong signal. 

Computation of the Code Dependent Crosscorrelation Factor 

The code dependent portion of the crosscorrelation factor is computed from 
the known relative states of the PRN code generators to predict the crosscorrelation 
between a strong signal of unit power and zero frequency offset, and a weak signal. 
This factor is multiplied by the amplitude of the corresponding strong signal and 
adjusted for frequency attenuation before it is subtracted from the composite signal. 

The various Gold codes used to modulate the PRN signals are all derived 
from a 2 code sequence G1 and G2 where the bits of the two sequences are 
combined through an XOR operation after G2 has been offset some number of bits 
relative to.G1 depending on the Gold code selected. As noted elsewhere in the 
specification, an XOR operation using binary numbers is mathematically equivalent 
to multiplication of ± 1. This allows expressing the equations below in term of 
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products of ± 1 while in reality the implementation could be with binary numbers with 
XORs, 

The correlation between two C/A codes can generally be expressed as 
ESat1G1(l) X Sat1G2(l) x Sat2G1(l-offset) x Sat2G2(l-offset) x e'^^®' 

where 

I = Summation index ranges from 0 to 1022 

Sat1G1(l) = Value of satellite Vs G1 coder chip at state I. Possible values are ±1 
Sat1G2(l) = Value of satellite Vs G2 coder chip at state I. Possible values are ±1 
Sat2G1(l) = Value of satellite 2's G1 coder chip at state I. Possible values are ±1 
Sat2G2(l) = Value of satellite 2's G2 coder chip at state I. Possible values are ±1 
offset = time difference between the satellite 1 and 2 in units of chips 
AG = Phase change per chip between satellite 1 and 2 in radians 

It should be noted that when the difference l-offset is less than 0, 1023 is . 
added to the difference to maintain the value in the range of 0 to 1022. In other 
words, the domain of the functions returning coder chip states is limited to the range 
of Oto 1022. 

The computation time required to compute the 1023 bit-by-bit correlations 
can be accelerated by making use of standard CPU instructions that perform 8, 16 
or 32 bit-wise XORs with a single CPU instruction. The following will demonstrate 
the method of computing 8 chips in parallel. Those skilled in the art will immediately 
recognize that the scheme can be easily modified to accommodate some other 
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convenient number of bits per CPU XOR operation. 

The 1023 states of G1 and G2 are stored linearly in permanent memory. 
Thus it is possible to quickly gather 8, 16, 32 or some other convenient number of 
bits with a single CPU load instruction by computing the address of the desired chip 
and the shift required to align it. Thirty-two bits is a particularly convenient number 
because 31 divides 1023 evenly. The preferred embodiment thus reads 32 bits at 
a time and uses 31 of them at a time for each of 33 intervals that span the 1023 
chips of the C/A code. The 31 bit sums are broken into 4 parts of 8, 8. 8, and 7 bits, 
and each 7 or 8 bit sum is multiplied by e"'^^' where I changes by 7.75 chips for each 
part. The form of the sum is 

2(g-Mei«3i ^ 2(Sat1G1(|x31+J) X Sat1G2(|x31+J) X 

Sat2G1(|x31+J-offset) x Sat2G2(|x31+J-offset) ) + 
Q-jAeo«3i*7.75)^ S(Sat1G1(|x31+J+8) x Sat1G2(|x31+J+8) x 

Sat2G1(|x31+J+8-offset) X Sat2G2(|x31+J+8-offset) ) + 
e-jAe(i«3i*i5.5)^ E(Sat1G1(|x31+J+16) x Sat1G2(|x31+J+16) x 

Sat2G1(|x31+J+16-offset) x Sat2G2(|x31+J+16-offset) ) + 
e-'^^<"'^'*"^5»x2(Sat1G1(|x31+J+24) x Sat1G2(|x31+J+24) x 

Sat2G1(|x31+J+24-offset) x Sat2G2(|x31+J+24-offset) ) ) 

where 

I = Outer index ranges from 0 to 32 

J = • . Inner index ranges from 0 to 7 for the first three sums- and from 0 to 6 for 
the last sum. The inner sums are computed in parallel by using a 32 bit 
word that contains all 31 bits and using bitwise XOR to perform the 
multiplications and shifting and adding to sum the 1 bit products. 
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Note that all of the multiplications of the G1 and G2 codes in the above 
equation are implemented by bit-wise XOR instructions. The above algorithm is in 
error by at most -17dB from an exact computation, and requires about 6000 CPU 
operations to complete. 

Use Of Computed Crosscorrelations 

Periodically, as needed, the code dependent crosscorrelation factors are 
computed for all strong and weak signal pairs with small frequency differences, i.e. 
frequency differences that could cause strong-weak crosscorrelation interference. 
In the preferred embodiment strong signals are those with C/Nq > 40 dB and weak 
signals are those with C/Nq < 30 dB. Because 10 ms integrations of I, Q 
measurements are used by the code and phase tracking loops, the maximum 
"significant" frequency difference (modulo 1000 Hz) is 90 Hz. In the preferred 
embodiment the code dependent cross correlation factor for each possibly 
interfering pair of signals is computed for each of the measurements that might 
potentially be used by the tracking and signal processing algorithms. For example, 
if eariy, punctual and late measurements are used by the tracking loops, the 
correlation factors for each of these code alignments is computed and stored in the 
tables. 

These tables only need be updated at a 10 Hz rate because the maximum 
Doppler difference is less than 9 KHz or less than 6 chips per second. In addition 
to maintaining the crosscorrelation table, the frequency attenuation of the 
crosscorrelation due to the frequency difference is computed at the 10 Hz rate. The 
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attenuation can be expressed as 



Frequency Attenuation = sin(AF mod 1000 x n/100)/(AF mod 1000 x n/100) 
where 

AF = Frequency difference between a strong and weak signal in Hz 

Mod = modulo offset to give a range of -500 Hz to +500 Hz 

The attenuation only needs to be recomputed if the frequency difference changes 
by more than 5 Hz. 

Strono Signal Estimation and Removal 

An estimate of the phase and amplitude of the strong signal is required to 
remove the crosscorrelation. The method used in the preferred embodiment is to 
track the strong signal on its own dedicated channel and collect the I. Q 
measurements output over the exact same Interval that the weak signal I, Q 
samples are taken. The known phase and frequency of the replica signal that is 
used to track the strong signal is an excellent approximation of the actual phase and 
frequency of the strong signal. Furthermore, because the strong signal is in phase 
lock, the magnitude of the I measurement provides a good approximation of the 
amplitude of the strong signal. Finally, the bl-phase modulation of the strong signal 
data bits D can cause the phase of the strong signal to- rotate 180 degrees 
whenever the data bits transition from 1 to 0 or from 0 to 1. In the preferred 
embodimenti the phase of the strong signal is corrected by adding 180 degrees to 
the phase of the replica signal whenever the sign of the I measurement for the 
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strong signal is negative. 

Every 10 ms a new set of I, Q correlation data is available from the channel 
assigned to track the weak signal. The tables of crosscorrelation factors are 
checked to predict the presence of any interfering strong signals. If strong signals 
are predicted, the following subtraction is performed to remove the strong signal 
crosscorrelations 



FirstCodeOffset = WeakCodeState - StrongCodeState - StrongDoppler xaT 

- TableEntryOCodeState 
DeltaPhase = WeakCarrierPhase - StrongCarrierPhase - StrongDoppler x AT 

+ DeltaKHz x StrongCodeState 
FirstPhase = FirstCorrelationPhase + DeltaPhase 
SecondPhase = SecondCorrelationPhase + DeltaPhase 
FirstMag = FirstCorrelationMag + FirstCodeOffsetFraction x StrongI 

X FrequencyAttentuation 
SecondMag = SecondCorrelationMag x (1 -FirstCodeOffsetFraction) x StrongI 

xFrequencyAttenuation 
Corrected WeaklQ = WeaklQ - FirstMag x g-i^'^s^P'^^s® - SecondMag x Q-iSecon^Phase 



where 

WeakCodeState = 
StrongCodeState = 
StrongDoppler = 
AT = 



Code state of last output to the weak signals channel 
Code state of last output to the strong signals channel 
Doppler of last output to the strong signals channel 
Difference in time between outputs to the weak and 
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strong channels 

TableEntryOCodeState = Code state difference of the first element of the 

crosscorreiation table. 
WeakCarrierPhase = Carrier phase angle of last output to the weak signal 

channel 

StrongCarrierPhase =Carrier phase angle of last output to the strong signal 

channel 

DeltaKH2= Nearest integer multiple of IKH2 of the difference 



of KHz. 

FirstCorrelationPhase = Phase entry in the crosscorreiation table for the chip 

indicated by FirstCodeOffset 



chip indicated by FirstCodeOffset + 1 chip. 
FirstCorrelationMag = Magnitude entry in the crosscorreiation table for the chip 

indicated by FirstCodeOffset. 
SecondCorrelationMag = Magnitudeentry in the crosscorreiation table for the chip 

indicated by FirstCodeOffset + 1 chip. 
FirstCodeOffsetFraction = Fraction of a chip in FirstCodeOffset. 
StrongI = Absolute value of I correlation from the strong channel. 

FrequencyAttenuation = Attenuation due to frequency offset. 
WeaklQ = |Q correlation from the weak signal's channel. 

CorrectedWeaklQ = IQ correlation corrected for crosscorreiation from the 



between the weak and strong channels Doppler. In units 



SecondCorrelationPhase 



Phase entry in the crosscorreiation table for the 



strong signal. 
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CorrectedWeaklQ is computed for the early, on time, and late correlators by 
shifting the FirstCodeOffset appropriately, such as by half a chip each. These 
modified correlations are then used normally in the carrier and code tracking 
software for the weak signal. The algorithm attenuates the crosscorrelation by at 
least 10 dB without attenuating the weak signal, and is repeated for each strong 
signal that may be interfering with the weak signal. 

While a particular form of the invention has been illustrated and described, 
it will be apparent that various modifications can be made without departing from the 
spirit and scope of the invention. It must be understood that although the preferred 
embodiment has been described in term of application to GPS receivers, the 
method of the invention can be applied to any other communication system 
employing CDMA spread spectrum transmissions that are susceptible to the near- 
far problem and where the interfering strong signals can be measured or estimated . 
with sufficient accuracy. Accordingly, it is not intended that the invention be limited 
except as by the appended claims. 



-24- 



Patent Application 
Docket No. SiRF-49352 



